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Abstract 

Within the massive field theoretical renormalization group approach the ex- 
pressions for the (3- and 7- functions of the anisotropic ran-vector model are 
obtained for general space dimension d in three-loop approximation. Resum- 
ming corresponding asymptotic series, critical exponents for the case of the 
weakly diluted quenched Ising model (m = 1, n = 0), as well as estimates 
for the marginal order parameter component number m c of the weakly di- 
luted quenched m-vector model are calculated as functions of d in the region 
2 < d < 4. Conclusions concerning the effectiveness of different resummation 
techniques are drawn. 
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Introduction 



Study of the critical behaviour of the Ising model has several attractions. On the 
one hand, the Ising-like models are simple enough, which is of a special advantage in 
the statistical physics. On the other hand, in spite of their simplicity, such models 
show rich and interesting behaviour at the critical point. Also, the existence of 
the exact solution for the two-dimensional Ising model often makes it an object for 
verifying different approximation schemes. All the stated above yielded the high 
interest devoted to the problem. In particular a great deal of generalization of the 
model appeared. Among different ways of generalization, much attention has been 
devoted to the affect of the impurities on the critical behaviour of the Ising-like 
models as well as to the investigation of critical regimes of the models on the lattices 
of a non-integer dimension (d). There have been devised different realizations of the 
last stated generalization. For example, one can approach the concept of non-integer 
dimensionality either by explicit construction of the non-integer dimensional object, 
which leads to the concept of a fractal [|J, or by formal carrying out an analytic 
continuation of the function, which by definition depends on a natural value of 
dimension. 

Within the theory of critical phenomena the latter ambiguity was reflected in 
examining the critical behaviour of the many-particle systems on fractal ||,[| or on 
abstract hypercubic lattices of the non-integer d. There arosed a question whether 
a model on a fractal lattice (being scale invariant) possesses universality as well as 
a system on a hypercubic lattice (having translation invariance). The problem has 
been widely studied but still remains open Q-Q. Today's point of view states that 
the usual demand for "strong universality" (in sense of critical properties depending 
only on symmetry of the order parameter, interaction range and space dimension) 
seems not to be obeyed by fractal lattice systems, and for them the concept of 
universality itself should be revised 

Speaking about the study of Ising-like models on analytically continued hypercu- 
bic lattices of non-integer d, one should note a great variety of theoretical approaches 
devised for these problems. These include: the Wilson-Fisher e-expansion |jlO| im- 
proved by the summation method pif ; Kadanoff lower-bound renormalization ap- 
plied to some special non-integer dimensions [[12j|; high-temperature expansion im- 
proved by a variation technique [13|; finite-size scaling method applied to numerical 
transfer- matrices data [14,1{|; new perturbation theory based on the physical branch 
of the solution of the renormalization group equation [16-19]; fixed dimension renor- 
malization group technique [20,21] applied directly to arbitrary non- integer d ppl. 
Perhaps the first paper devoted to the study of the Ising model in different, however 
not non-integer dimension, was [24] where non-universal properties of the model were 
discussed. All these approaches, as well as the computer simulations, confirm the 
correctness of the universality hypothesis also for non-integer d hypercubic lattices 
and allow to obtain the critical exponents as functions of d with high accuracy. 

In spite of the variety of approaches to treat the general non-integer d case 
[11-19, 22, 23] the results for the Ising model critical exponents obtained on their 
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basis lay close to each other; the above mentioned analytic continuations may appear 
actually equivalent. Note, however, that the analytic continuation in e = 4 — d at 
general non-integer dimension leads to the fail of the Yang-Lee theorem ]25| . On the 
other hand, the study of Ising-like spin systems on non-integer dimension hypercubic 
lattices cannot be reduced to the fractal lattice systems [25] except for the case of 
vanishing lacunarity limit H and thus is the task of individual interest. 

Returning to the study of the critical behaviour at integer d, one should note 
that the problem becomes more complicated when studying spin systems with a 
structural disorder. Whereas the case of the annealed disorder is of less interest 
from the point of view of determining asymptotical values of critical exponents [^] , 
the weak quenched disorder has been a subject of intensive study. Here the Harris 
criterion |2^| has been devised. It states that if the heat capacity exponent a pure of 
a pure model is negative, that is the heat capacity has no divergence at the critical 
point, impurities do not affect the critical behaviour of the model in the sense that 
critical exponents remain unchanged under dilution. Only in the case a pure > 0, 
the critical behaviour of the disordered model is governed by a new set of critical 
exponents. As far as for a 3d m- vector spin model only the 3d Ising model {m = 1) 
is characterized by a pure > 0, it is the Ising model which is of special interest. And 
because of the triviality of the annealed disorder in the sense mentioned above, the 
most interesting object for study is just the quenched Ising model. The appearance 
of a set of new critical exponents for that model at d = 3 is confirmed by the 
experiments 2£-31|, renormalization group (RG) calculations [32-41], Monte-Carlo 
(MC) |42|-[46| and MCRG (4| simulations. 

The situation is not so simple for the 2d Ising model. Onsager exact solution 
of the pure model proves the logarithmic divergence of heat capacity, which yields 
oipure = 0, and allows one, in accordance with the Harris criterion, to clasify this 
case as a marginal one. Most of the theoretical works suggest that the 2d Ising 
model with a quenched disorder has the same critical behaviour as the 2d pure Ising 
model (except for logarithmic corrections) []39|,^C|,^-|53|] (see also review |54j]). This 
result is corroborated by MC-simulations on two-dimensional lattices [55|-|59[ and 



exp er iments [30,61]. 

Deviations from the expected critical exponents, which sometimes are observed 
during such computations, are explained by a system being not in the asymptotic 
region (see E^] for recent study). Nevertheless, some authors assert that for the 2d 
Ising model with a quenched disorder a new critical behaviour appears 62, 3j| 

While the undiluted Ising model at non-integer d was a subject of intensive 
study 



10 19], it is not the case for the diluted Ising model. Only the work []36| can 
be mentioned here, where the model was studied within the Golner-Riedel scaling 
field 1 64] approach. It is worthwhile to note that the e-expansion technique applied 
to this model, due to the fact that .RG-equations appear to be degenerated on the 
one loop level, results in -y/e-expansion for the critical exponents [34|. The latter 
is known up to the three-loop order |3^, ^6| . The equations of the massive field 
theory at fixed integer d [2C, 21] first applied to the diluted Ising model at d = 2,3 
in [|35,[37 were found to be the most effective method for investigating this problem. 
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In order to consider an arbitrary non-integer d the Parisi approach [EH, 21] was 



generalized in [67] where critical behaviour of the model was studied in a two-loop 
approximation. The aim of the present work, based on the massive field theoretical 
approach, is to make a more detailed investigation of the critical behaviour of the 
diluted 0(m)-vector model at arbitrary d. Though it is the case to = 1 in which we 
are interested most of all, we consider the -RG-equations for any m, which also allow 
us to study the crossover in the model at any d. We will obtain the -RG-equations 
within the 3-loop approximation and apply to their analysis different resummation 
procedures in order to find the most reliable one. 

The set-up of the article is as follows. In the next Section we introduce the 
model and the notation. Then we describe the i?G-procedure adopted here and give 
the series for the i?G-functions of the weakly diluted quenched m-vector model in 
the three-loop approximation. Being asymptotic, these series are to be resummed. 
This is done in Section 2 where different ways of resummation are used. Section 3 
concludes our study giving results for the quantitative characteristics of the critical 
behaviour and discussing them. In the Conclusions we give some general comments 
to the present work. In the Appendix we list some lengthy expressions for the 
coefficients of the i?G-functions in the three-loop approximation. 



1 The Model and the RG - procedure 

As it is well known, the critical behaviour of the quenched weakly-diluted m - vector 
model is governed by a Lagrangian with two coupling constants [54]: 
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in replica limit n — > 0. Here any <f> a is a m-component vector 

(j) a = ((ft 01 ' 1 , (f) a > 2 , . . . ,c/> a ' m ); no > 0,vo < are bare coupling constants; too is bare 
mass. 

As it was already stated above, we adopt here the massive field theory renor- 
malization scheme [20|, 21] in order to extract the critical behaviour governed by 
([!]) . We start from the defined by (|l|) unrenormalized one-particle irreducible vertex 
functions 

r (L,Ar) (gi, .., q L ;pi, ..,p N ; to , u , v ; A ; d) (2) 

depending on the wave vectors {q}, {p}, bare parameters too, uq, vq and the ultravio- 
let momentum cutoff Ao- The vertex functions' dependence on the space dimension 
d is explicitly noted here as well. We impose the renormalization conditions at zero 
external momenta and non-zero mass (see [68, 6S] for instance) at the limit Ao — ► oo 
for the renormalized functions [7C] r^' 2 \ r^^^T^'^ , r^' 2 ^' 
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r i?' 2) (P> -p;m,u,v;d)\ p=0 = m 2 , (3) 
^2 r ^' 2) ^' -P;m,u,v;d)\ p=0 = 1, (4) 

r i?,« ({P*}; m ^ u i d ) l{Pi}=o = m A ~ d u, (5) 
rg; 4) (te};m,u, V ;d)| {pi}=0 = m^v, (6) 

r fl'^(g , ;p ) -p;m-, , ",w;d)| 9= p = o = 1, (7) 

with m,u,v being the renormalized mass m = Z 3 mi = Z 3 r(°' 2 )(0; mo,uo,vo) and 
couplings u = m d ~^Z^Zi u uo, v = m d ~^Z 2 Z^v§. From these conditions there follow 
expansions for the renormalized constants for field (Z3), vertices u (Z\ u ), v {Z\ )V ) 
and cj) 2 insertion {Z?). Subsequently, these define the coefficients /?,7 entering the 
corresponding Callan-Symanzik equation: 

du 

Pu{u,v) = — \ UQV0 , (8) 

a mm 

dv 
a mm 

7</> = 73 = \uo,v , (10) 

a mm 

In the stable fixed point {u*,v*} to be defined by simultaneous zero of both 
/3-functions: 

u (u*,v*) = 0, 

(3 v (u*,v*) = 0, (12) 
the 7^-function gives the critical exponent 77 of the pair correlation function: 

l4> ( U *,v*)= V . (13) 
The correlation length critical exponent v is defined in the stable fixed point by: 

^{u*,v*) = 2-v- x - l4> {u*,v*). (14) 

Using familiar scaling relations, one can easily calculate any other critical expo- 
nents on the base of r\ and v. 

Applying the described above procedure, one obtains in the three-loop approxi- 



mation [[n] (3- and 7- functions in the form [72\ 
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Here d is the space dimension, to is the order parameter component number, n is 
the replica index, i\ and 72 are dimensionally dependent two-loop integrals. The 
corresponding coefficients for three-loop parts are listed in the Appendix. The val- 
ues for the three-loop integrals 73 ... is which appear in three-loop coefficients for 
integer d = 2, 3 are listed in |fT3[ . In particular, substituting loop integrals 71,72 
as well as 73, ... ,78 in (JT5|) - (|l~8|) by their values at d = 3 we get at n = 0,m = 1 
the corresponding functions of the 3d weakly diluted Ising model, which in the 3- 
loop approximation were obtained in |35|]. At d = 3, to, n- arbitrary corresponding 
expressions coincide with those, obtained for the 3d anisotropic mn-vector model 
in [74|. Our idea is to keep the dimensional dependence of the loop integrals and, 
being based on their numerical values for arbitrary d |23||, to study the 0(mn)-model 
at arbitrary (non-integer) d as well. But for the reason explained above, the point 
of main interest here will be the replica limit n = of the anisotropic mn-vector 
model, especially the case m = 1. 

Expressions for (5- and 7-functions will be the starting point for the qualitative 
study of the main features of the critical behaviour which will be done in the next 
section. 



2 The Resummation 



As we have already mentioned, the values of the 7-functions in a fixed point (u*,v*) 
lead to the values of the critical exponents 7/ and v. However, it is well known now 
that the series for KG- functions are of asymptotic nature [75-77] and imply the 
corresponding resummation procedure to extract reliable data on their basis. Let us 
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note, however, that as to our knowledge the asymptotic nature of the series for RG- 



functions have been proved only for the case of the model with one coupling [78], and 
the application of a resummation procedure to the case of several coupling constants 
is based rather on general belief than on a proved fact. One of the resummation 
procedures, which in different modifications is most commonly used in the studies 



of asymptotic series, is known as the integral Borel transformation [79|. However, 
this technique implies explicit knowledge of the general term of a series and thus 
cannot be applied here, where only truncated sums of the series are known. To get 
over this obstacle one represents the so-called Borel-Leroy image of the initial sum 
in the form of a rational approximant and in such a way reconstitutes the general 
term of the series. The technique which involves a rational approximation and the 
Borel transformation together, is known as the Pade-Borel resummation technique 



(in the field-theoretical RG content see (81,82j as an example of its application). 



Note here that the resummation technique, based on the conformal mapping, 
which is widely used in the theory of critical phenomena [pj|, cannot be applied in 
our case because its application postulates information on the high order behaviour 
of the series for (3- and 7-functions. The latter is still unknown for the theory with 
the Lagrangian ([[]). 

To summarize up the stated let us write that the Pade-Borel resummation is 
performed as follows: 

• constructing the Borel-Leroy image of the initial sum S of n terms: 

S = Ya i x i ^Y °' (at)t - (19) 

where T(x) is the Euler's gamma function and p is an arbitrary non-negative 
number. The special cases p = and p = 1 correspond to resumming (5- 
functions without or with prefactors u and v in accordance with the structure 
of the functions (I!)-©; 

• the Borel-Leroy image ( |i~9| ) is extrapolated by a rational approximant 
[M/N] (xt), where by [M/N] one means the quotient of two polynomials; M 
is the order of the numerator and N is that of the denominator; 

• the resummed function S res is obtained in the form: 

S res = / dtexp(-t)t p [M/N](xt). (20) 
Jo 

In the two variables case only the first step is changed; namely, here we define 
the Borel-Leroy image as 

0<i+j<n 0<i+j<n 1 ^ ^ J ^f^r 1) 

Generalization to the many-variable case is trivial. 
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Now one can take into account that the second step of the stated scheme can be 
done in different ways. One can write down various Pade approximants in the vari- 
able t to obtain within the three-loop approximation the expressions of the structure 
[2/1], [1/2] and [0/3]. It is also possible to use Chisholm approximants [84] in the 
variables u and v, which, generally speaking, in the same number of loops can be 
of type [3/1], [2/2], [1/3] and [0/4], but the explicit definition of any approximant 
needs some additional equations now P4|. The technique, which involves Chisholm 
approximation together with the integral Borel transformation is referred to as the 
Chisholm-Borel resummation technique. To be consistent, one would have to apply 
the all different resummation frameworks in order to obtain reliable results on their 
basis and find which of the methods is the most effective. However, strong restriction 
on the number of choices can be imposed. 

First of all, an approximant should be chosen in the form reconstituting the sign- 
alternating high-order behaviour of the general term of (3- and 7-functions, which 
was confirmed in the particular case m = 1, n = 2 and n = 3 [pq| . The approximant 
generating a sign-alternating series might be chosen in a form [M/l] with the positive 
coefficients at the variable t (or u and v ). Choosing an approximant with a non-linear 
denominator, generally speaking, one does not ensure the desired properties. Direct 
calculations affirm the argumentation: /?- functions, resummed with the Pade-Borel 
and the Chisholm-Borel methods with approximants [M/N] ,N>1, for u < 0, v > 
give the roots which lie far from the expected values which for d = 3 are known up 
to the order of four loops ]39| and for general d were calculated from the two-loop 
/3-functions [S7|. This is true for any p. The stated results permit us to eliminate 
from the consideration approximants with a non-linear denominator. 

Note as well that choosing representation of the i?G-functions (|i~5|)-(|i~8|) in the 
form of Pade or Chisholm approximant of type [M/l] might result in the appearance 
of a pole in the expression. Here we use an analytical continuation of the resulting 
expressions by evaluating the principal value of the integral. Treating the task in this 
way one notes that the topological structure of the lines of zeros for the resummed by 
the Pade-Borel technique /3-functions is very different in the region near the solution 
for the mixed fixed point and strongly irregular when passing through the number 
of loops. In particular this yields that in the three-loop approximation there exist 
two solutions close to the expected value of the mixed fixed point. To compare, the 
results obtained within the frames of the Chisholm-Borel method do not have these 
faults and are more stable from the point of view of proceeding in number of loops. 

So, the results given below are obtained by the Chisholm-Borel method applied 
to the approximant of type [3/1]. In order to determine the form of this approx- 
imant completely one must define two additional conditions. The approximants 
are expected to be symmetric in variables u and v , otherwise the properties of the 
symmetry related to these variables would depend, except for the properties of the 
Lagrangian, on the method of calculation. By the substitution v = all the equa- 
tions which describe the critical behaviour of the diluted model are converted into 
appropriate equations of the pure model. However, if pure model is solved indepen- 
dently, the resummation technique with the application of Pade approximant is used. 
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Thus, Chisholm approximant is to be chosen in such a way that, by putting any of 
u or v equal to zero, one obtains Pade approximant for a one-variable case. This 
also implies a special choice of additional conditions. In the present study amidst 
all the possible expressions which satisfy the stated demand we choose Chisholm 
approximant [3/1] by putting coefficients at v? and i> 3 to be equal to zero. 



3 Results 

Now we are going to apply the mathematical framework which was discussed in 
previous sections in order to obtain numerical characteristics of the critical behaviour 
of the weakly-diluted Ising model in general dimensions. It was noted in the Section 1 
that the critical behaviour of the quenched weakly-diluted Ising model is described 
by the effective Lagrangian (|l|) in the case m = 1 and zero replica limit n = 0. 
Namely, the task in the end comes to obtaining fixed points which are defined by 
simultaneous zero of the both /3-functions. Among all the possible fixed points one is 
interested only in those in the ranges u* > 0, v * < and only in stable ones where the 
stability means that two eigenvalues 61,62 of the stability matrix B = d(3 Ui /ditj\ u * , 
Ui = {u, v} are positive or possess positive real parts. The structure of the /3- 
functions ji~5|)-(|l6|) yields the possibility of four solutions for the fixed points. The 
first two {u* = 0, v* = 0} and {u* = 0, v* > 0} in our case at d < 4 are out of 
physical interest, while the second pair which consists of pure {u* > 0, v* = 0} and 
mixed {u* > 0,v* < 0} points, are responsible for two possible critical regimes. 
The critical behaviour of the diluted model coincides with that of the pure model 
when the pure fixed point appears to be stable. If the mixed point is stable, the 
new (diluted) critical behaviour of the system takes place. The type of the critical 
behaviour depends on the number m of the order parameter components and on the 
dimensionality d: at any d, 2 < d < 4 a system with large enough m is not sensitive 
to the weak dilution in the sense that asymptotic values of critical exponents do 
not change; only starting from some marginal value m c , at m < m c a mixed fixed 
point becomes stable and the crossover to the random critical behaviour occurs. The 
problem of determining m c as a function of d will be discussed later. Now we would 
like to state that m c > 1 for any d, 2 < d < 4, and thus just the mixed fixed point 
governs the asymptotic critical behaviour of the diluted Ising model. 

If one attempts to find the fixed points from the /3-functions ([i"5|)-(ll6|) with- 
out resummation, there always appears only the Gaussian {u* = 0, v * = 0} trivial 
solution; the existence of the rest possible three fixed points depends on the con- 
crete details of the /3-functions portions in the braces in expressions (|D^)-(|l6|). In 
a 3d case it appears that without a resummation the non-trivial mixed fixed point 
does not exist in one-, two- and four- loop approximations [E9lEo|. It is only the 



three- loop approximation where all the four solutions of the set of equations (12) 
exist [35]. In figure [l] we show the behaviour of the non-resummed /^-functions of 
the three-dimensional weakly diluted Ising model in the three-loop approximation. 
Resummed functions are shown in the same approximation in figure H. Note that in 
this approximation the shape of the functions remains alike in the region of small 
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couplings u and v. Fixed points correspond to the crossing of the lines f3 u = 0, /3 V = 
as it is demonstrated in figures ||[ |I[ The left-hand column in figures H |I] shows the 
lines of zeros of non-resummed /3-functions in three-dimensions in one-, two-, three- 
and four-loop (results of |4(|) approximations. One can see in the figures that 
without resummation all non-trivial solutions are obtained only within the three- 
loop level of the perturbation theory. In the next order all fixed points disappear 
which is a strong evidence of their accidental origin. At any arbitrary d, 2 < d < 4 
the qualitative behaviour of the functions is very similar to that shown in figures |3| 
and |. 

As it has already been mentioned, in order to reestablish the lost pure and 
mixed points one applies the resummation procedure to /3-functions. In the three- 
dimensional space the result of resummation is illustrated by the right-hand column 
in figures ||] and |||. Here we have used the Chisholm-Borel resummation technique 
choosing Chisholm approximant in the form discussed in the previous Section with 
p = 1 in successive approximation in the number of loops. The icons in the figures 
which correspond to a one-loop level are the visual proof of the degeneracy of the 
/3-functions in this order of the perturbation theory: the plots of root-lines are paral- 
lel independently of resummation. The rest three images in the right-hand columns 
are a good graphic demonstration of the reliability of the Chisholm-Borel resum- 
mational method: two-, three- and four-loop pictures are quantitatively similar, the 
coordinates of the pure and mixed point are close. 

The numerical results of our study are given in table |]. Here, the coordinates of 
the stable mixed fixed point and the values of the critical exponents of the quenched 
weakly diluted Ising model are listed as functions of d between d = 2 and d = 3.8. 
The eigenvalues b\ and 62 of the stability matrix are given as well. 

It was already noted that the values of 7-functions in a stable point yield the 
numerical characteristics of the critical behaviour of the model. For example, given 
the resummed functions 7^ es and 7^2 es , the pair of equations 

7 f = V, (22) 

T$>{u\v*) = 2-u- l - V (23) 

allows us to find the exponents r] and v. All other exponents can be obtained from 
the familiar scaling laws. 

However, one can proceed in a different way. That is, by means of the scaling 
laws it is possible to reconstitute the expansion in coupling constants of any expo- 
nent of interest or of any combination of exponents, and only after that to apply 
the resummation procedure. If exact calculation were performed the answer would 
not depend on the sequence of operations. However, this is not the case for the 
present approximate calculations. We have chosen the scheme of computing where 
the resummation procedure was applied to the combination zy 1 — 1 = 1 — 7^2 — 7^ 
and 7 _1 = (2 — 7^2 — 7^)/(2 — 7^). The exponents a, j3 and rj have been calculated 
on the basis of numerical values of the exponents 7 and v. The resummation scheme 
appears to be quite insensitive to the choice of the parameter p given by fll9|), fl20|). 
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However note, that computations have been performed here, as well as in [67|, with 
p = 1. 

Comparing our data from table |l] for the critical exponents at d = 2 with the 
results for the pure Ising model one can see that the exponent 7 differs from the 
exact value 7/4 by the order of 5%, the exponent v is smaller from the exact value 
v = 1 less than by 4%. This confirm the conjecture that the critical behaviour of 
the weakly diluted quenched Ising model at d = 2 within logarithmic correction 
coincide with that of the pure model (see [54] for review). It is also interesting to 
compare numbers given in table |l] with those obtained for general d within the 2-loop 
approximation |67j : all the exponents of the three- loop level lie slightly farther from 
the expected exact values of Onsager than those of the two-loop approximation. 
This may be explained by the oscillatory nature of approaching to the exact values 
depending on the order of the perturbation theory. It is also interesting to note 
that the two-loop approximation yields better estimates for the heat capacity critical 
exponent a for all d in the range under consideration. Namely, in accordance with the 
Harris criterion, the exponent a for the diluted Ising system should remain negative. 
This picture is confirmed much better by the two-loop approximation where a is 
negative in the whole range of d, unlike the three-loop level of the perturbation 
theory, the results of which yield a > for 2 < d < 2.8. 

However, table |l] shows that the next (third) order does improve our understand- 
ing of the critical behaviour of the model in general dimensions. The results of the 
two- loop calculations [S7| show that starting from some marginal space dimension 
the approach to the stable point becomes oscillatory: the eigenvalues 61 and 62 turn 
to be complex possessing positive real parts. This is an artifact of the calculation 
scheme and therefore it was expected |37|] that by increasing the accuracy of calcu- 
lations one decreases the region of d which corresponds to the complex eigenvalues. 
It is really the case. In the three- loop approximation the region of complex 61,62 
is bounded from below by d = 3.3, whereas in the two-loop approximation |)7j the 
corresponding value is lower and is equal to d = 2.9. Thus, the region of d character- 
ized by the oscillatory approach to the stable fixed point shrinks with the increase 
of the order of the perturbation theory. 

The comparison of the three-dimensional value of v with the four-loop result p(J 
v = 0.6701 gives the accuracy of 0.05% for our computations (compare with 1% for 
two-loops, where the value u{d = 3) = 0.678 was obtained). Thus, it may be stated 
that the general accuracy of calculations decreases when passing from d = 4 to d = 2 
which, in particular, results from the fact that our approach is asymptotically exact 
at upper critical dimension d = 4. 

The comparison of the present results with the other data available is provided by 
figure |5[ Here, the behaviour of the correlation length critical exponent v obtained 
by different methods is demonstrated in general dimensions. The results of the 
massive field-theoretical scheme are plotted by solid (three-loop approximation; the 
present paper) and dashed (two- loop approximation; ref. (67|) lines. One can see 
that the two lines practically coincide far enough from d = 2, in particular, both 
lie very close to the most accurate result for d = 3 [EJ which is shown by the box. 
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The application of the scaling-field method [36| yields numbers shown in figure |5] by 
stars. The limit from below (d = 2.8) of the method applicability is caused by the 
truncation of the set of scaling-field equations, which was considered in [|36|| . One can 
also attempt to obtain some results by resumming the y^-expansion which is known 



for the diluted Ising model up to three- loop order [65, S6 and for the exponents v 
and 77 reads: 

1 1/6 U/2 535 - 756C(3) 
" = 2 + 4(o3 £ ) + 8(53)* £ ' (24) 

* = -4 + (l) 1/2 w( 24+7 ^ 3 )^ 3/2 ' ™ 

where ("(3) ~ 1.202 is Rieman's zeta function. The corresponding results are shown 
by open diamonds. They were obtained by applying the Pade-Borel resummation 
scheme to the series of y/e- expansion (]24]) |65| , |66| . The value of v obtained in 
such a way is of physical interest only very close to d = 4. Even in the next 
orders of the expansions the values of critical exponents are not improved |^q |; 
this is an evidence of the -y/e-expansion unreliability in tasks like the one under 
consideration. To compare, one can state that the situation with the applied in the 
present paper theoretical scheme is contrary to the y^-expansion. While the two- 
loop approximation is valid in ranges 2 < d < 3.4, the next order of the perturbation 
theory enlarges the upper bound up to d = 3.8. One can expect that the next steps 
within the perturbation theory will allow one to obtain the description of the critical 
behaviour of the model with enough accuracy for any d, 2 < d < 4. 

Let us recall now that expressions (]T5|)-(|l8|) for the i?G-functions, as well as 
their three-loop parts listed in the Appendix, allow us to study asymptotic critical 
properties of the mn- vector model with arbitrary m and n in arbitrary d not only 
for the case m = 1, n = 0. In particular, by keeping m as an arbitrary number 
and putting n = one can obtain the numerical estimates for the marginal or- 
der parameter component number m c which divides the diluted (governed by the 
mixed fixed point) asymptotic critical behaviour from the pure one, when the 0(m)- 
symmetric fixed point remains stable. In accordance with the Harris criterion the 
case m = m c corresponds to zero of the heat capacity critical exponent a of the 
model. One may extract the value of m c from this condition. However, the above 
discussed results of the three-loop approximation do not yield enough accuracy for 
a. Alternatively, the fixed mixed point should coincide with the pure fixed point at 
m = m c , which in particular means that v*(m = m c )\ m i xe d = 0. The last condition 
was chosen as a basis of our calculation. The appropriate numbers of the present 
three-loop approximation (thick solid line) together with the data of the two-loop 
approximation (dashed line) |67| are shown in figure |6|. The result of e-expansion 
m c = 4 — 4e is depicted by the thin solid line. In the three-loop approximation we 
obtain m c = 1.40, d = 2 and m c = 2.12, d = 3. These values are to be compared 
with the exact results of Onsager which yield m c = 1 at d = 2, and the theoretical 



estimate m c = 1.945 ± 0.002 [37|. One can see that the two-loop results are closer 



to the expected values for both d = 2 and d = 3. For a two-dimensional case the 
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two-loop value m c = 1.19 [67] differs from the exact one by 20%, while the three- 
loop number decreases the accuracy to 40%. The case d = 3, m c > 2 contradicts the 
suggestion that the xy-model asymptotic critical behaviour should not change un- 
der dilution in three-dimensions. The reason for decreasing the calculation accuracy 
with increasing the order of the perturbation theory may lie in oscillatory approach 
to the exact result. One can expect that already the four-loop case will improve the 
estimates for m c for all 2 < d < 4. Let us also note that the determination of m c 
may serve as a test for improving the resummation scheme. 



4 Conclusions 

The goal of this paper is to study the critical behaviour of the weakly diluted 
quenched Ising model in the case when the space dimension d continuously changes 
from d = 2 to d = 4. 

As it was mentioned in the Introduction, the study of the pure Ising model at 
arbitrary d, which corresponds to a scalar field-theoretical model with one coupling 
constant, is the subject of a great deal of papers. It is not the case for the model 
with a more complicated symmetry. In particular, here we study a model with 
two couplings corresponding to terms of different symmetry in the Lagrangian ([!]). 
Such a problem was studied previously on the basis of the scaling-field method p6|], 
and field-theoretical fixed dimension renormalization group calculations within a 



two-loop level of the perturbation theory are available [S7|. 

Our calculations hold within the theoretical scheme of [23, 37 . This approach 
appears to be one amidst other possible calculation schemes for many tasks; however, 
in our case it seems to have no alternatives within the field-theoretical approach. 

Being asymptotic, the resulting series for the i?G-functions are to be resummed. 
In the present study we have chosen the Pade-Borel and the Chisholm-Borel resum- 
mation techniques. Restricting ourselves to analytic expressions for the resummed 
functions, we present numerical data mainly obtained on the basis of the Chisholm- 
Borel resummation technique. Note that the absence of any information on the 
high-order behaviour of the obtained series for the -RG-functions does not allow one 
to apply other resummation schemes, e.g. those based on the conformal mapping 
technique |j83[ . 

The quantitative description of the critical behaviour of the model is steady from 
the point of view of passing from the two- to the three-loop approximation. Smaller 
agreement between the two- and the three-loop approximations at d far away from 
d = 4 may be explained in a way that the precision of computing falls down with 
the increase of the expansion parameter which takes place at decrease of d. The real 
parts of eigenvalues corresponding to the mixed point seem to remain positive up to 
d = 4, which testifies that at arbitrary d the weakly diluted quenched Ising model 
is described by the mixed fixed point. 

The work was supported in part by the Ukrainian Foundation of Fundamental 
Studies (grant No 24/173). 
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Appendix 



Here we have collected the most lengthy expressions for the three-loop contributions 
to the i?G-functions. The three-loop part of the /3 u -function reads: 

df A (u, v) = ffi°u 3 + Pl x u 2 v + /ft V + /ft V, (26) 

where 



Pt 



■3,0 



Pu 



■2.1 



PL 



3° 



.0,3 



4(31m 2 + 430m + 1240)n + (m + 8)(m + 2) x 



(m + 8) 3 

( - (3d + 8)z 2 + 12(i 3 + i 8 )) + 48(m 2 + 20m + 60)i 4 + 
24(2m 2 + 21m + 58)i 5 + 6(3m 2 + 22m + 56)i 6 + 
24(5m + 22)i 7 + 8(4m 2 + 61m + 178)1 ; 

12(17m 2 + 256m + 780)n + 



1,2 _ 



(m + 8) 2 (mn + 8) 
(m + 2) ( - (3dm + 42d + 16m + 80)i 2 + 12(m + 14)i 3 + 

18(m + 8)i 8 ) + 24(3m 2 + 70m + 224)i 4 + 

6(15m 2 + 158m + 448)i 5 + 6(3m 2 + 32m + 100)i 6 + 

48(5m + 22)i 7 + 6(9m 2 + 146m + 448)1 ; 

12(19m 2 n + 80mn + 470m + 2032)ii 



(m + 8)(mra + 8) 2 
(s{mn + 8) (3d + 4) + m(3dmn + 40mn + 78d + 176)) i 2 + 

12(m 2 n + 8mn + 26m + 64)i 3 + 48(m 2 n + 8mn + 68m + 
292)i 4 + 12(llm 2 n + 34mn + 136m + 584)i 5 + 
6(m 2 n + 8mn + 50m + 256)i 6 + 576(m + 5)i 7 + 
36(m + 2)(mn + 8)i 8 + 12(5m 2 n + 22mn + 136m + 584) 

—, ^— T \e(mn + 10)(-2(mn + 23)n + 3i 6 ) + 

(mn + 8) 6 L 

(mn + 2) ( - (4mn + 9d + 8)z 2 + 36i 3 + 3(mn + 8)i 8 ) 
72(3mn + 22)i 4 + 9(m 2 n 2 + 14mn + 88)z 5 + 
24(mn + 14)z 7 + 3(m 2 n 2 + 38mn + 264) . 

The three- loop part of the /3„-function reads: 

fi™(u, „) = ^3^3 + 01,2^2 + p 2,l u 2 v + / g3 > u 3 j 

where 



+ 



(27) 



18 



ft 



0,3 



1 



(mn + 8) 3 
(mn + 2 



- 4(31mV + 430mn + 1240)ii + (mn + 8) 



ft 



1,2 



ft 



2,1 _ 



ft 



3,0 



■ (3d + 8)i 2 + 12(t 3 + *s)) + 48(m 2 n 2 + 

20mn + 60)i 4 + 24(2m 2 n 2 + 21mn + 58)i 5 + 6(3m 2 n 2 + 

22mn + 56)i 6 + 24(5mn + 22)i 7 + 8(4m 2 n 2 + 61mn + 178) 

4(m + 2) r 

-7 ,7/ r - 4(28mn + 275)n - 

(mn + 8) 2 (m + 8) L v Ji 

(3dmn + 4mn + 15d + 56)i2 + 12 (mn + 5)i 3 + 

24(2mn + 27)i 4 + 3(13mn + 100)z 5 + 6(3mn + 13)i 6 + 

96i 7 + 9(mn + 8)i 8 + (29mn + 316)1 ; 

-7 — 77— ^ \ ~ 12(mn + 42m + 224)ii - 

(mn + 8)(m + 8)^ L 

(3dmn + 12dm — 8mn + 48d + 16m + 256)^2 + 
12(mn + 4m + 16)z 3 + 48(5m + 34)i 4 + 
12(13m + 56)i 5 + 6(3mn + 14m + 40)i 6 + 
144i 7 + 36(m + 8)i 8 + 12(llm + 64) 
2(m + 2) r 



m + ■ 



4(llm + 70)n - 3(dm + 2d + 16)i 2 + 

6(m + 2)(2i 3 + 3i 6 ) + 2(m + 8)(12i 4 + 3i 5 + 3i 8 + 5) 
The three-loop part of the 7^-function reads: 



lf A {u,v) 



m + 2 
(m + 8) 2 
3(m + 2 



u° + 



3(m + 2) 



(m + 8)(mn + 8) 

9 mn + 2 o 
uv + — V 



u 2 v + 



(m + 8)(mn + 8) (mn + 8) 2 

The three-loop part of the 7U2 -function reads: 



(3t 8 - 4*2). 



, \ -3,0 3 1 -2,1 2 , -1,2 2 , -0,3 3 

b 2 (u, v) =7^ti + y 2 u v + 7^ uv + y 2 v , 



where 



-3,0 



-2,1 

7^,2 



- 4(llm + 70)n + (m + 2) ( - (3d - 8)i 2 + 



m + 2 

(m + 8) 3 

12i 3 + 18i 6 ) + 2(m + 8)(12i 4 + 3i 5 + 5) 

- 12(mn + 10m + 70)n + 



m + 2 



(m + 8) 2 (mn + 8) 
(mn + 2m + 6)( - (3d - 8)i 2 + 12i 3 + 18i 6 ) + 
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6(ra + 8)(12i 4 + 3i 5 + 5) 
3(m + 2) 



4(llmn + 70)n + 
(m + 8)(mn + 8) 2 L 

(mn + 2)( - (3d - 8)z 2 + 12t 3 + 18i 6 ) + 

2(mra + 8)(12i 4 + 3« 5 + 5)l; 

mn+2 r „. . 

— - - 4(llmn + 70)n + 

(mn + 8y L 

(mn + 2)( - (3d - 8)i 2 + 12i 3 + 18i 6 ) + 
2(mn + 8)(12i 4 + 3i 5 + 5) 
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FIGURE CAPTIONS. 
Figure The non-resummed /3-functions in the three-loop approximation; d = 
3, m = 1, n = 0. The dark surface corresponds to the /3 u -function. 

Figure Q The Chisholm-Borel resummed /3-functions in the three-loop approxima- 
tion; d = 3,m = l,n = 0. The dark surface corresponds to the /3 u -function. 

Figure || The lines of zeros of non-resummed (left-hand column) and resummed 
by the Chisholm-Borel method (right-hand column) (3- functions for m = 1, n = 
in different orders of the perturbation theory: one- and two-loop approximations. 
Circles correspond to (3 U = 0, thick lines depict V = 0. Thin solid and dashed 
lines show the roots of the analytically continued functions (5 U and f3 v respectively. 
One can see the appearance of the mixed fixed point u > 0, v < in the two-loop 
approximation for the resummed /3-functions. 

Figure |I| The lines of zeros of non-resummed (left-hand column) and resummed by 
the Chisholm-Borel method (right-hand column) /3-functions for m = 1, n = in 
three- and four-loop approximations. The notations are the same as in figure |3| 
Close to the mixed fixed point the behaviour of the resummed functions remains 
alike with the increase of the order of approximation. This is not the case for non- 
resummed functions. 

Figure ||. The correlation length critical exponent v of the weakly diluted Ising 
model as a function of the space dimension d. The results of two- |37j] and three- 
loop (the present paper) approximations are shown by the dashed and the solid lines 
respectively, the square reflects the number of the four-loop approximation [40| at 



d = 3, stars correspond to work [|36]] and open diamonds refer to the resummed 
-\/i-expansion. 

Figure |6[ The dependence of the marginal order parameter component number m c 
on the space dimension d. Two- and three-loop results are shown by the dashed and 
thick solid lines respectively, the e-expansion data m c = 4 — 4e are depicted by the 
thin solid line. 
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Figure 6: 

Table 1: The stable point coordinates, critical exponents and the eigenvalues of 
the stability matrix of the weakly diluted Ising model at arbitrary d. The three- 
loop approximation (the superscript " c" denotes that real parts of the corresponding 
eigenvalues are given). 
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